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SUMMARY 



This investigation was performed to determine the feasibility of 
obtaining lateral frequency response zind stability derivatives utilizing 
the forced oscillation technique. Tests were conducted in a North 
American Navion at one altitude, airspeed and gross weight condition. 

Both rudder and ailerons were oscillated sinusoidally with a 
simple mechanical oscillator and aircradEt responses were telemetered 
and recorded on magnetic tape. A technique was developed for elimin- 
ating distortion of the inpvit sine wave by utilizing a low frequency band 
pass filter in the reduction of flight data. Frequency responses were 
determined for (p , ~\p' , and ^ for both the rudder and the 

aileron oscillations; and these were used to determine the aircraft's 
stability derivatives. Frequency responses were calculated from the 
experimentally determined derivatives and were compared with the 
actual frequency responses to illustrate the accviracy of the method 
employed to extract the stability derivatives. 

The methods employed proved entirely feasible for this aircr<ift 
at the single flight condition and apparently were applicable to many 
current production aircraft. The frequency responses showed little 
raindom error, eind the responses calculated from the experimentally 
determined derivatives closely matched the actual responses. 
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INTRODUCTION 



The design of modern high speed aircraft and the automatic 
control systems with which nearly all are equipped requires an 
acciirate evaluation of their response characteristics. In the early 
design stages the only method available to the designer involves the 
estimation of the eiircraft's stability derivatives and the solution of 
the assiimed equations of motion. Because of the inherent difficulties 
involved in the analytic determination of stability derivatives, large 
errors may be present in the derived response characteristics. 
Dynamic flight testing of the prototype aiircraft to accurately determine 
its stability derivatives and the coefficients of its transfer fxmction is 
becoming increasingly important. 

Dynamic flight test data is usually obtained by one of two 
methods. The most popular of these involves the determination of the 
transient response of the aircraft to known control inpiits. While the 
data thus obtained is, in itself, not in the best form for analysis, the 
recent development of efficient computing inachines has made it possi- 
ble to convert such a transient response to frequency response form 
where it may be more readily analyzed. The second method involves 
the direct determination of the aircraft's frequency response by the 
steady state forced oscillation technique. While considerably more 
flight time is required in this method, data reduction is greatly simpli- 



2 



fied; and good res\alts may be obtained without the use of automatic 
computers. 

The determination of an eiircraft's lateral frequency response by 
the forced oscillation technique may be accomplished by oscillating 
either of the lateral controls with a sinusoidal forcing function while 

keeping the other locked. From simultaneous recordings of the air- 

* * 

craft responses in (p , 1/^ and ^ and the control inputs on 

separate time traces, it is possible to determine amplitude ratios eind 
phase angles for each response over a range of oscillation freqxaencies 
for the steady state condition. From these amplitude and phase spectra 
it is, in theory, possible to determine from either the aileron or rudder 
data all the 2 iircraft*s aerodynamic derivatives, except of coxirse the 
control derivatives for the locked surface, providing the weight and 
moments of inertia of the aircraft are known. In practice, however, it 
is frequently difficult to determine all the derivatives accurately since 
some of them have little effect on the frequency response and thus are 
difficult to determine from it. 

While considerable effort has been expended on this type of testing, 
especially at the Cornell Aeronautical Laboratory, the acctirate deter- 
mination of derivatives without harmonic analysis of the time traces has 
required very pvire sinusoidal control inputs. This requirement has led 
to complicated instrumentation which has made the method applicable 
only to research type testing. 
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This investigation was conducted to determine the applicability 
of the method to more simplified testing techniques utilizing a simple 
mechanical oscillator, telemetering eqviipment for data recording, and 
an tiltra-low frequency band -pass filter to eliminate sine wave distor- 
tion. The flight tests were conducted in a North American Navion at 
the Forrestal Research Center, Princeton University, during the 
winter and spring of 1957. 
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EQUIPMENT 



Test Airplane 

The test airplane was a North American Navion, which is an all 
metal, four place, low wing, single engine airplane. The Navion is 
equipped with a Continental E-185 engine, rated at 185 HP at 2300 
RPM and 29 inches of MP, at sea level. The wing is xmswept, and the 
tedl configuration conventional. The control surfaces are ciU of beaded 
skin construction, eind have no trailing edge extensions. Trim tabs 
are located on the elevator, edleron, and rudder. 

The test airplane was modified by the removal of the two rear 
cockpit seats to provide space for the test instrumentation. The test 
airplane had a gross weight of 2870 pounds with fxill service, pilot and 
co-pilot, and all test equipment. 

A three view drawing of the Navion is shown in Fig. 1 and the 
general specifications of the test airplane are tabulated in Table I. 

Fig. 2 is a photograph of the test airplane. 

Power Supply 

Power available for the instrvimentation in the test airplane in- 
cluded both 12 and 28 volt direct current supplied by a 28 volt engine 
driven generator and two 12 volt batteries connected in series. A 500 
VA voltage and frequency controlled inverter was installed behind the 
co-pilot's seat to provide 115 volt, 400 cycle, three phase alternating 
ctirrent power to the gyros and telemeter transmitter. 
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Considerable effort was expended in obtaining an inverter which 
could be frequency controlled, since it was very important that the 
three phase power to the rate gyros be constant at the calibration 
frequency of 400 cycles. 

A circviit diagram of the aircraft 28 volt system and the inverter 
system is shown in Fig. 3. 

Aircraft Instrumentation 

Standard aircraft instruments were used to determine indicated 
airspeed and altitude, hieasurement of the following quantities was re- 
qxiired: sideslip angle, yaw rate, and roll rate with respect to wind or 

stability axes, and aileron and rudder position angles. In choosing the 
type of components needed to measure the above quantities, an attempt 
was made to utilize instrumentation compatible with the telemetering 
equipment. For this reason, the transducers selected were of the 
potentiometer type, not only because they were capable of providing 
signals of sufficient magnitude, but also because their excitation source 
readily provided zero and full scale reference voltages. 

^ view of some past history of various amounts of coiq)ling be- 
tween the response of the aircraft and certain of the instrumentation 
components in this type of flight testing, considerable attention was 
given to choosing transducers with satisfactory frequency response 
characteristics in order to avoid correction of response data for the 
dynamics of the recording instruments. In general, it was desirable 
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that the natural frequencies of the rate gyros and the sideslip vane be 
large as compared to the expected forced oscillation frequencies of the 
aircraft. 

Sideslip angle was sensed by a light vane moxinted on a boom 
approximately one maximvim fuselage diameter ahead of the left wing 
tip. The moment of inertia of the vane was very low, and no gearing 
was introduced in the linkage in order to keep the natural frequency and 
damping within acceptable limits. An Electro-Mechanical Laboratory 
ultra-low torque potentiometer with .008 inch-ounces of friction torque 
was obtained for use with the sideslip vane. 

In order to reduce boom resonance, which had been experienced 
in previous tests with the Navion, the boom was stiffened by three angle 
strips equally spaced anoimd the periphery, and running the full length 
of the boom. The sideslip vane and the boom are shown in the assem- 
bled condition in Fig. 4. 

Aileron and rudder angle were measured by 10, 000 ohm Model 
5301 single tvirn Helipot precision potentiometers. Due to the nature of 
the differential aileron system of the Navion, it was not considered de- 
sireable to pick off aileron angle at the aileron itself. Accordingly, the 
pickoffs were made from the control cable movement just inboard of 
the right aileron. A position close to the aileron was chosen in order 
to avoid the effects of cable stretch. The rudder deflection potentiome- 
ter was mo\onted on the bottom of the rudder and was actuated by a 
cable under spring tension. 
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The measviretnent of yaw and roll rates was accomplished with 
type 361Z Giannini rate gyros. This type of gyro operates on a power 
supply of 115 volt, 400 cycle, three phase AC to the motor, and 115 volt 
single phase AC to the heaters. The transducer is of the potentiometer 
type, with a DC resistance of 5250 ohms. These gyros were obtained 
on loan from the Instriunentation Laboratory of the U. S,^ Naval Air 
Test Center, NAS Patuxent River, Maryland. 

Oscillating Mechanism 

Design of the control driving mechanism was based on the follow- 
ing criteria: (a) The system must generate a sinusoidal output of con- 
trollable frequency which may be applied to either the rudder or the 
aileron control; (b) The driving link for the driven control, and the 
locking mechanism for the undriven control must each include a quick 
disconnect, not subject to binding when under load; (c) Driving and lock- 
ing links must include a trim adjustment; (d) Inter -unit rigidity and 
close tolerances are required due to the cyclical varying load; and (e) 
The unit must be compatible with the power and space limitations of the 
Navion. 

These conditions were met in the following manner. The rotation 
of a 400 cycle, l/2 HP DC motor was converted to a pure sinusoidal 
translational motion by a scotch yoke, with suitable intermediate gear- 
ing. The driving rod to the rudder could be disconnected from the 
scotch yoke plate and bolted to the mounting plate to provide a positive 
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rudder lock. The rod from the scotch yoke plate was fitted inside a 
tube from the rudder. By use of a large handled set screw the rod and 
tube could be rigidly locked at any desired trim position, and quickly 
vmlocked without possibility of binding. The aileron drive was obtained 
by leading light aircraft control cable around pvilleys from the horizon- 
tally translating shaft of the scotch yoke to a set of vertically translat- 
ing turnbuckles, from which a chain drove a sprocket wheel moiinted 
on the control coltuxm. The sprocket was secured to the control colvimn 
by a large handled set screw arrangement which allowed aileron trim 
adjustment and a quick disconnect. To lock the ciilerons, the drive was 
opened at the turnbuckle and a cable bolted to the base plate substituted 
for the driving cable. 

The motor, gear box, and driving and locking unit was built up on 
one large section of 1/4" aluminum plate, which provided a very rigid 
base. The complete mechanism could be installed or removed as a unit, 
the plate being bolted to extrusions riveted to the siffeners \inder the 
deck on the starboard side of the cockpit. In-flight shifting from the 
aileron drive to the rudder drive involved abovit five minutes. 

Control of the driving motor was obtained by the use of an ampli- 
dyne movinted behind the co-pilot. Through control of the field of the 
amplidyne, whose output provided the armature current to the driving 
motor, any speed from zero to maximum rated, including reversing, 
could be obtained. Indication of the operating frequency was not directly 
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available in the cockpit, but could be obtained approximately from a 
calibrated voltmeter across the amplidyne signal circuit. At the lower 
frequencies, the torque (proportional to the armature current) became 
marginal, and was evidenced as a distortion to the sine wave of the 
driven control as the aerodyiiamic load varied with control displacement, 
A worm and gear was incorporated in the gear box to obtain irreversi- 
bility between the motor and the control svurfaces, but the variation of 
the frictional load on the worm was sufficient to cause this low frequency 
distortion. The aircraft instrumentation is shown in Fig. 6 and a close 
up of the oscillating mechanism is shown in Fig. 7. 

Signal Circuit Design 

The telemeter transmitter was designed to accept voltages from 
the transducers from zero to five volts. It was thus desirable for the 
full scale voltages on the potentiometers to be maintained at five volts. 
Variations in this full scale value were permissable, providing that the 
variation on all the potentiometers was the same, since the telemeter 
corrected all voltages to a percentage of the full scade reference voltage. 
The potentiometer used in the sideslip vane was a 360 degree single 
turn model, and, since only 40 degrees of variation in sideslip angle were 
expected, it was necessary to apply 45 volts excitation to this potentiome- 
ter in order that its oxitput might have a five volt variation. The use of 
a 45 volt DC signal circuit supply was thus dictated. It was also desira- 
ble to maintain the signal voltage to the other potentiometers at five volts 
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and introduce a biasing circuit to provide a variation of sideslip output 
voltage from zero to five volts rather than a five volt variation about 
some intermediate voltage. 

The signal circtiit design based on the preceding requirements is 
shown in Fig. 5. The essential element of the circ\iit is the bias sup- 
ply voltage potentiometer, pot ”A". With the sideslip vajie against the 
stop in one direction, pot "A" was adjusted so that the output voltage 
of the sideslip pot was zero. The vame was then deflected to the oppo- 
site stop, and the output voltage of the sideslip pot was recorded. 

This value was slightly less than five volts. Pot "B" was then adjusted 
to bring the ftill scade reference voltage, and thus also the voltage in 
the remaining signal circviits in agreement with the fiall scale sideslip 
voltage. The procedure was repeated two or three times as necessary, 
in order to nullify any loading on the system. Once the final adjustment 
was made, the zero to full scale voltage for each quantity to be meas- 
ured was exactly the same, and small variations in signal voltage pro- 
dviced the same variation for all circvdts. Pots "A" and "B" were 
locked and required no further adjustment for the remainder of the 
project. 

IXuring the initial flight and transmission tests, considerable AC 
ripple was fo\xnd to enter the DC sign 2 il voltages, very likely from the 
400 cycle inverter. Condensers "C* and *D" in Fig. 5 were installed, 
and were found to reduce the ripple to a negligible value. 
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Telemeter System 

The telemeter unit was of the pulse width type and was designed 
and built by the Applied Science Corporation of Princeton (ASCOP). 

This unit was the key element in the instrxunentation of this project. 

The basic functional components of this pulse width system are trans- 
ducers (all of the potentiometer type); a time division multiplexer, or 
commutator; an electronic cod.er, or keyer to put the intelligence into 
piilse form; an RF transmitter; and groimd station eqmpment for re- 
ceiving, decoding, and recording the data in usable form. 

The airborne component of the system is designated the ASCOP 
D Series PW Multicoder, and the ground station equipment the ASCOP 
M Series PW Ground Station. 

A unique feature and very great advantage of pulse width data 
systems is that, due to the mxiltiplexing or time sharing operation, all 
data channels traverse the same path during transmission. Since zero 
and full scale references are transmitted with the data, this results in 
a calibration signal, effective for all channels, being transmitted 20 
times each second. These calibration signals allow the removal of 
almost all of the system error, with the exception of errors introduced 
in the instrument pickups and the mechanical switch, and they make 
possible the use of this telemetering system in tests of long diiration 
smd varying environmental conditions. 
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Magnetic Tape Recorder 

An Ampex Model 309 C dtial track three speed tape recorder was 
an integral component of the ground station equipment, and provided a 
convenient means for recording, storage, and playback of the standard 
ptilse width modulated data. Provisions were also included for placing 
voice transmissions directly on the magnetic tape from either or both 
the aircraft and the ground station. 

Data Recorder 

Qualitative analysis of the flight data during the actual flights and 
qTiantitative analysis of the recorded signals was made on a Sanborn 
Model 154-lOOD four channel recorder. This instrument had a very 
high natural frequency (42 cps), and the low frequency response was 
flat to zero cps. Different paper speeds and variable gain were possible 
and permitted expansion of the time trace and amplification of the signal 
to nearly einy desired value. 

Filter 

Li order to study the response of the aircraft to a true sinusoidal 
driving function, all forcing functions and the corresponding responses 
were filtered to remove all but the fvindamental frequency. The filter 
employed was the model 330-A KROHN-HITE Ultra-Low Frequency 
Band-Pass Filter. This is an adjtistable filter whose gain is unity in 
the pass band and drops oiitside the pass band at the rate of 24 db per 
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octave. The filter combines a high pass filter and a low pass filter, 
each of which may be continuously tuned from 0.02 to 2, 000 cps. 

Since the filter will handle an input signal of only up to 10 volts rms 
without overloading, a voltage divider was applied to the telemeter out- 
put, which had a full scale range of up to 100 volts. For this applica- 
tion the filter was used as a narrow pass band with the high and low 
dials set at the same frequency. Used in this manner, the filter did 
not alter the relative amplitude or phase relationships between the 
input and the outpvit at a specific frequency. 

A photograph of the ground instrumentation is shown in Fig. 8. 

Instrviment Calibration 

Six items required calibration in the instrumentation system. 
These were the two rate gyros, the two control deflection potentiome- 
ters, the sideslip vane potentiometer, and the signal voltage to the 
amplidyne which was roughly a measure of the frequency of control 
oscillation. The aircraft altitude and airspeed indicating systems were 
not calibrated because previous calibrations indicated that the errors 
here were negligible for the conditions of this investigation. 

The two rate gyros were calibrated on a pendulum which consisted 
of a 3/4“ diameter aliiminum shaft and a 12" diameter instrviment pan. 
The pendulum, the overall length of which was 121.75" was pinned to 
the overhead truss structure in the Forrestal hangar. A scale was 
mounted with brackets on a small table in such a way that displacement 
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of the pan from its rest position could be measured. The gyros were 
mounted one at a time on a wooden base plate which was clamped to the 
pendulvim pan, and their potentiometer outputs were connected to a 
Brush Recorder. The pendulvun was displaced from its rest position 
and allowed to swing freely. Since there was an excessive amount of 
friction in the system, swing amplitudes were measured on two succes- 
sive swings in the same direction and averaged. The two maximum 
voltages recorded during the swing were also averaged and the maximum 
rate attained during the swing was determined from the following 
relation: 

A _ d 3(^0 

“max p 

where P was the natural period of the pendulum whose measured val^xe 
from the trace on the Brush recorder was 3. 33 sec. , 1 was the length 
of the pend\ilxim, and d was the average maximum displacement of the 
pan from its rest position. The pendulum was allowed to swing until 
once again at rest with readings taken, as described above, over the 
entire range of gyro rates. The angular rate gradient of both gyros 
was found to be linear and equal to . 610 degrees/sec. per percent full 
scale potentiometer voltage. The gyro calibration curves are shown in 
Fig. 9. It was found that the gyros were very sensitive to the frequency 
of the 400 cycle supply voltage and as a resxolt great care was taken 
during calibration amd also during the actual flight test program to 
maintain the frequency of the supply at exactly 400 cycles per second. 
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The control deflection potentiometers were calibrated by measur- 
ing the resistance between the output terminals over the full range of 
control deflections. The actvial control deflections were measiired with 
a vernier protractor capable of accuracy to 0.1 degree. The rudder de- 
flection gradient was found to be linear and equcil to .405 degrees per 
percent fiiU scale potentiometer voltage. The Navion aircraft is 
equipped with differential ailerons and it was expected that the aileron 
Cedibration would therefore be non-linear. However, using average 
aileron angle, it was found that over the small deflection range of +10 
to -10 degrees desired during the flight test program the edleron deflec- 
tion gradient was linear and equal to .1933 degrees per percent full 
scale potentiometer voltage. 

The sideslip vane potentiometer was calibrated in the same man- 
ner as the control deflection potentiometers, using a protractor and 
measuring the resistance between the outpvit terminals with an ohm- 
meter. The calibration curve is shown in Fig. 10. The sideslip 
gradient was foxind to be linear and eqtocil to .468 degrees per percent 
full sccile potentiometer voltage. 

The telemeter equipment corrected its output to percent fviU scale 
voltage and since full scale voltage on all the potentiometers varied by 
the same increment as the reference voltage for small fluctuations of 
the latter value, it was xmnecessary to calibrate the potentiometers 
through the telemeter eqviipment. 

The speed control on the oscillating device was calibrated roughly 
in flight by timing the frequency of control oscillations and recording 
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the voltmeter reading across the control field of the amplidyne genera- 
tor. This speed control was used only to obtain order of magnitude 
settings, the actual frequency of any particvdar oscillation being deter- 
mined from the trace on the Sanborn Recorder. 

Careful consideration was given to the possible errors which 
might be introduced in the data by the dynamics of the recording instru- 
ments. The natural frequency of the rate gyros was 17 cycles per 
second and their damping ratio was between .45 and .7. The natural 
frequency and damping ratio for the sideslip vane were experimentally 
determined in flight to be 15.44 and 0.081 respectively. Since the 
natural frequency of the aircraft was approximately 0.4 cycles per 
second and the maximum forced oscillation frequency was not expected 
to exceed one cycle per second, the attenuation of the signals due to the 
dynamics of the instruments was effectively zero. While a small error 
in phase angle measvirement was introduced at the higher frequencies, 
it was considerably smaller them the reading accviracy with which phase 
angle cotild be determined from the data, and thus was not considered 
in the reduction of flight test data. 
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FLIGHT TEST PROCEDURE 



In this investigation, frequency response was measvtred at the 
single flight condition of 5000 ft. density altitude, 129.6 mph true air- 
speed, cind 2870 lbs. gross weight. Each test was commenced at 
steady level flight with the amplidyne control set at zero speed, the 
oscillating mechanism exactly centered to its mid position, and the 
lateral controls locked in this trim position. The run was then initi- 
ated by setting the speed control and allowing steady state conditions 
to become established. The necessity for driving the lateral controls 
about the exact trim position was so great that the lock-on system to 
a moving oscillator, as employed in a somewhat similar longitudinal 
study (Ref. 1) was not feasible. In spite of the great care exercised 
to center the oscillations about a straight and level flight path, prac- 
tically every run was viltimately terminated by the excitation of the 
spiral mode. At this time the rapid unlocking features in the me chan - 
ical system design were most appreciated. 

In setting the speed control, ein approximation to the desired 
frequency could be obtained by the use of a voltmeter across the ampli- 
dyne control field calibrated to read frequency. The relatively short 
d\iration of each run did not permit fine adjustments to obtain a specific 
frequency, so several runs around a desired frequency (primarily near 
the dutch roll) were employed. In addition, it was found that the low 
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frequency runs were best started by going rapidly to a high frequency, 
and then coming down to the desired frequency. Otherwise, the 
response to the slow initial deflection drove the aircraft so far into the 
spired divergence that the subsequent control reversal was insufficient 
to establish an oscillation about a straight flight path. Although it had 
been expected that cdl controls except the driven one would be locked 
d\iring a run, it was found that slight elevator movements were required 
to maintain a constant airspeed throughout the run, particularly dxiring 
oscillations of the rudder . 
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RESULTS 



Data Reduction 

The qviantities measiired d\iring the flight test program were Sa 
• • 

ox t fS > <p t and . These four quantities were teleme- 

tered to the ground station, recorded on magnetic tape, and simulta- 
neoxisly plotted on the Sanborn Recorder. This raw data had noise auid 
distortion superimposed on the fxindamental simxsoidal driving frequency. 
An example of such a record is shown in Fig. 11. In order to study the 
response of the aircraft to only the fundamental of the driving frequency, 
a KROHN-HITE Ultra -Low Frequency Band-Pass Filter was used to 
filter both the driving function and the aircraft responses. From the 
original record of the raw data the fundamental frequency for a particu- 
lar run could be determined with sufficient accuracy to establish the 
high and low settings for the band pass filter. With the filter thus set 
to pass only the specific frequency, all four quantities were successively 
played from the magnetic tape through the telemeter receiver (convert- 
ing the pvilse width coding to a varying voltage output), thence through 
the filter to the Sanborn recorder. As each qxiantity was being filtered 
and recorded, an unfiltered trace of the forcing fvinction was also 
placed on this new record to act as a time reference. The results of 
this time and recorder paper consuming process were four sinusoidal 
records of the measiired quantities, each with a common unfiltered time 
reference. Fig. 12 shows an example of such a filtered record with the 
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unfiltered forcing functions matched for each response. From these 
matched responses the amplitude ratios and phase angles relative to 
the driving qxxantity could be determined. Since the settings on the 
filter and the recorder were not altered during the filtering of a given 
frequency, any amplitude attentviation or phase shifting was applied 
equally to all foxir quantities. The gain of the recorder was initially 
set for each frequency to accommodate the quantity with the greatest 
amplitude, such that this quantity would utilize almost full scale on a 
recorder paper channel. By using this technique, and increasing the 
paper speed for the higher frequencies, a rather accurate determina- 
tion of the amplitude ratios and phase angles was possible. 

The quantities 0 and were determined from the records 

* a 

of <p and by integration. This was readily performed by divid- 

ing by u) and subtracting 90° from the phase angle, since the re- 
sponse was sinusoidal. The magnitude of each response was divided 
by the inagnitude of the forcing function and the appropriate scale fac- 
tors were applied to give the desired amplitude ratios. Since the sideslip 
vane was not located at the c. g. of the aircraift zm additional correction 
to the /S’ response was required. This correction was a function of 
the distance of the sideslip vane from the c.g. of the aircraft and is 
given by the following relation; A . oZ 7 4- ’Y' 

The final frequency responses were tabulated in Tables II and in and 
were plotted in Figs. 13-18 and labeled "experimental." 



21 



Also plotted on the appropriate graphs in Figs. 13-18 were fre- 
quency response curves denoted as ’•theoretical." These curves were 
the result of an initial study of the Navion airplane based on lateral 
stability derivatives determined from previous tests, such as described 
in Ref. 2, or from theoretically determined derivatives. These initial 
"best guess" derivatives are tabiilated in Table IV, and were used as a 
rough check on the resvJts, as well as the basis for the initial system 
design. 

Deteirmination of Derivatives 

Ref. 3 develops a method for extracting stability derivatives from 
freqxiency response data for the lateral case. The method Involves the 
application of the vector technique to the three linear equations of lateral 
motion to separate them into six equations. The vector responses are 
resolved par£dlel and perpendicxilar to the control deflection vectors. 

The perpendicxilar resolution eliminates the control deflection deriva- 
tives in one of the two eqviations derived from each basic eq\iation and 
thus permits solution of the aerodynamic derivatives separately. The 
development of these equations is shown in Appendix A. The method of 
least squares is then applied to the data for each of these equations to 
yield lateral stability derivatives. 

Although the method of least squares is probably the best way to 
eliminate random error, it is not readily applicable to hand calculation 
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techniques. Automatic computing facilities were not available here, 
and since time was limited, it was decided to use the simpler method 
of averages in conjvinction with the vector equations developed in Ref. 3 
to determine the lateral stability derivatives of the Navion. The K 
notation of Ref. 3 was modified to a starred derivative notation for 
greater clarity. This notation is defined in Appendix A. The method 
of averages, as described in detail in Ref. 4, is also outlined in 
Appendix A. 

The components of the responses along and perpendiculaLr to the 
control deflection axis are defined in the form A+iB with the appropri- 
ate subscripts. Form\alas for the determination of the responses ajid 
their derivatives in this form are given in Appendix A. 

The two equations developed from the basic side force eqiiation 
by the vector technique are eqviations (13) and (14) in Appendix A. 

The equation describing the relation between the components perpen- 
dicular to the control deflection axis was: 



mine its value from the responses to either an aileron or rudder forcing 
function at a single frequency. Since the responses in yaw amd sideslip 
were greater for a rudder input, it was felt that the rudder data wo\ild 




Since the only unknown was it was possible, in theory, to deter- 
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the responses at 12 selected frequencies were taken from the faired 
rudder curves. The summation of the 12 eqviations determined from 
these 12 points gave one average equation, from which the value of 
Cy^* was determined. This gave a value of Cy^ of -0.591. The same 
technique using aileron data resulted in a positive value for , 
which was, of covurse, impossible. Ref. 3 mentions that the accurate 
determination of Cy^ may be difficult from this equation, since it 
involves the s\im of the yaw and sideslip components, which are nearly 
equal but opposite in sign. The axithors of Ref. 3 recommend obtaining 
lateral acceleration data and using a modified side force equation to 
avoid this condition. The large responses in yaw and sideslip to rtidder 
inpvcts tended to minimize this condition, but the determined value of 
Cy^ was nevertheless subject to considerable error. 

The equation describing the relation between the components 
parallel to the control deflection axis and therefore including the deriva- 



tive , C 



y/r 



was: 



(14) Cy^ - Ci_ Alp ~ +■ Sy.) 

Using the previously determined valoie for Cy^ in an average eqviation 
for the 12 selected points yielded a value for Cy^ This gave a 

value for Cy^^ of 0.206. The accuracy of Cy^^ was clearly depend- 
ent on the correctness of Cy^ , and was, in addition, svibject to small 
differences between yaw and sideslip components. Cy ^ was therefore 
of doubtful reliability. 
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The vector technique gave equations (15) and (16) of Appendix A 
from the basic roll equation. The relation between the components 
perpendicular to the control deflection axis was; 

(15) 4- 

Since was very small compared to Ix the term — was ass\imed 
to be zero. Solution for the three unknowns reqviired three equations. 
The 12 frequency points from the aileron data were now grouped to in- 
clude four frequencies below the dutch roll, four near the dutch roll, 
and four above the dutch roll, and average eqviations were determined 
for each group by adding the appropriate coefficients. Solution of these 
equations for Clp*, and were converted to , 

Cip, and Ci^, which were -0.0505, -0.436, and 0.0390, respectively. 
In the sol\ition of these equations it was foimd that an accvirate determin- 
ation of Ciy was difficxilt, since the solution involved a small differ- 
ence between large quantities; and sxziall errors were magnified many 
times. An attempt was made to carry greater accviracy in the solution 
than the data actually warranted, but of course this procedxure, while 
reqxiired to give any solution at all, cotild not give an accurate solvttion 
for Clp. This condition was not sxurprising, siitce it is well known 
that the cross derivatives are difficult to determine from this sort of 
data. 



Soltttion of the saune equation \ising a grouping of points from the 



pi lilt 
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rudder data should have given the same derivatives as previously 
determined from the aileron data. When this solution was attempted, 
however, it was foimd that the simultaneous equations were ill-condi- 
tioned. This arose from the fact that (9 emd 0 varied in nearly 
the same manner over the entire frequency range, both in phase angle 
and magnitude. Additional frequency points were taJcen and various 
groupings were tried in an attempt to eliminate this difficulty, but all 
were unsuccessful. Ref. 3 suggests a possible method of avoiding the 
diffictilty by eliminating the cross derivative coefficient by resolution 
of the vectors along the axis, instead of eilong and perpendic\ilar 
to the control deflection axis. This was attempted, but solution of the 
modified equation was still impossible. 

R was felt that the rudder data was slightly better them the aileron 
data, due to refinements in the flight test technique as the program pro- 
gressed, and it was therefore highly desirable to include the rudder 
data in the solution for the roll derivatives. The solvttion of the aileron 
eqtiations gave a value of Clp which was felt to be quite accurate, 
since it was in close agreement with theory and previoxis experimentsd 
results. This valxie was assvuned known in the rudder eqxiations for the 
rolling motion, and the two remaining derivatives Ci^ and Ci^ 
were determined by taking two groups of six frequency points each and 
arriving at two eqiiations in the two desired unknowns. The values de- 
termined for Cl- and Cx„ were -0. 0513 and 0. 119 respectively. 

P * 
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It was interesting to note that the value of Ci^ determined in this 
manner was nearly the same as that determined from the aileron data, 
and the doubtfvil derivative, Ciy, was nearer the theoretically deter - 
mined value. The derivatives determined in this manner were arbi- 
trarily chosen as the final values, since the authors felt that they repre- 
sented the best possible combination of all the experimental data. 

The relation between the components parallel to the control 
deflection axis was: 

(16) - Ctp u^^cp + yJ Ad) 

Using the previously determined values of and Cip* 

and utilizing one average equation based on all 12 frequency points 
from the rudder data permitted solution for which gave a 

of 0. 0128. Similarly, solution of the same equation using the 
aileron data gave Cl^ *, which resulted in a value for of 

-0.108. 

£bq>ressing the basic yaw equation in vector form resulted in 
eqviations (17) and (18) of Appendix A. For the components perpen- 
dicular to the control deflection axis the relation was: 

(17) ~ to ^ 

Because of the basic similarity between the equations for the yawing 
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ajid rolling motions, the vmknowns in equations (15) and (17) have the 
same numerical coefficients, with only the constant term differing. 
Because of this, rapid solution of the three average equations for the 
yawing motion was possible vising the same grouping of frequency 
points from the responses to aileron oscillations. This gave values of 

0.0771, -0.129, eind -0.0841, respectively. 
The same difficulty existed here in the determination of as had 

previously existed in the determination of Ci^. 

Because the numerical values of the coefficients of equation (17) 
were the same as those of equation (15), the three simultaneous equa- 
tions for the yawing motion derived from the rudder data were also ill- 
conditioned for solution. The method applied to solve equation (15) 
using an apparently accurate value of could not be applied in this 

case, because both auid as determined from the aileron 

data were felt to be doubtful, since they were not in agreement with 
either theory or the results contained in Ref. 5. In this case^final 
solution was effected by solution of a set of simultaneous equations, 
two of which were determined from the aileron data at frequencies near 
and above the dutch roll, and the third representing the average of all 
the rudder data. This approach appeared legitimate since each equa- 
tion in the set was independent of the type of forcing function. This 
particular choice of equations resulted from the opinion of the authors 
that the aileron data was better at higher frequencies and that inclusion 



28 



of the generally more reliable rudder data was desirable. It should be 
emphasized that this particular choice was peciiliar to this experiment. 
Solution of these eqxiations gave values for C 

of 0.0895, -0.106, and -0. 0908 respectively. 

The equation relating the components parallel to the control 
deflection axis for the yawdng motion was: 

(18) -Cy^^ + Cnp to Bet + ^- toB-f " ~ 

Using the previously determined values of 

and vitilizing all 12 frequency points from the rudder data permitted 
solxition for *, which gave a -0.0673. Similarly solu- 

tion of the same equation using the adleron data gave Cn *. which 
resulted in a value for of 0.00376. 

The final derivatives determined in this investigation are tabulated 
in Table IV. This table also contains the original theoretically deter- 
mined derivatives, and the derivatives determined for the same aircraft 
by an analog computer solution as outlined in Ref. 5. 

In order to check the accviracy of the method used to determine the 
aircrzift's stability derivatives, the transfer fxinction coefficients based 
on them were determined as OTJtlined in Appendix A, and from these the 
frequency responses were calculated. The experimental transfer fimc - 
tions are tabulated in Table V, and the frequency responses based on 
them are shown in Fig. 19-24. 
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DISCUSSION 



Dynamic flight testing is becoming increasingly importeint with the 
advent of high speed high altitude aircraft becatise the aircraft designer 
is more concerned than ever before with the accurate determination of 
stability derivatives and the associated transfer fimctions. Although 
transient methods are the most popular at the present time, it was felt 
that if the instrumentation reqviirements of steady state dynamic flight 
testing could be reduced eind the necessity for absoltttely sinusoidal con- 
trol inputs without the added requirement of harmozdc analysis could be 
eliminated the forced oscillation technique for determining frequency 
response and thus stability derivatives and transfer function coefficients 
could still be used to advantage xmder certain circvimstances. It was the 
purpose of this investigation to examine the possibility of obtaining ade- 
quate frequency response with relatively simple instrumentation, and if 
possible to eliminate the necessity of harmonic analysis due to distorted 
sinusoidal inputs. 

Since investigation of the methods involved was of primary con- 
cern here, flight tests were conducted at only one altitude, airspeed, 
and gross weight condition, although it should be emphasized that in any 
actucil flight test program to determine derivatives, all of these condi- 
tions should be varied over the operating range of the aircraft. 

Before the design of the instrumentation system could be accom- 
plished a prediction of the results based on theoretically determined 
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derivatives and control gradients was required. Obviotisly the experi- 
mental data had to include frequencies on either side of the dutch roll 
with frequencies at the dxitch roll being of primary interest. Previous 
investigations of the lateral dynamic stability of the Navion indicated 
that the period of the dutch roll was between 2 and 3 seconds. An 
analysis based on the best available derivatives, some calculated theo- 
retically from Ref. 6 and others taken from Ref. 2 and 7 was performed. 
This indicated that the period of the dutch roll was about 2.4 seconds. 
Theoretical frequency responses to both rudder and aileron forcing 
fxmctions were calculated and these served to design the instrumentation 
system. 

It was decided to design the oscillator to cover frequencies be- 
tween 0. 5 and 10 radians per second. The upper limit proved higher 
than necessary and better results would have been obtained if the design 
had limited the upper frequency to 5 radians per second. This upper 
limit was found to be well above the dxitch roll frequency and oscilla- 
tions at higher frequencies produced only negligible responses. Exam- 
ination of the control deflections necessary to produce recordable 
responses without exceeding the linear limit of the equations and the 
control gradients led to the decision to use a ^ H.P. motor and an 
amplidyne to control its armature current. Early experimentation with 
a mechanical Vciriable reduction gear indicated that this method was 
unsatisfactory. 
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Ref. 1 recommends a careful preliminary study of the dynamics 
of the recording instruments since, if their nat\xral frequency is not 
much larger than the natural frequency of the aircraft, large errors 
arise which must be corrected for. This was done, and it was found 
that the natural frequency of the rate gyros and sideslip vaine were such 
that the attenuation of their signal outputs was effectively zero and that 
the phase shift introduced was less than 0. 5°. Since other errors in 
phase angle measurement were expected to be of considerably greater 
magnitude, no attempt was made to correct for the dynamic response 
of the recording instruments. 

The instrumentation system based on the results of the above cal- 
culations proved completely adequate with the exception of the oscillator. 
It was found that since the output torque was directly proportional to the 
armature current of the motor the torque at low frequencies and thus 
low armature currents was insufficient to produce anything approaching 
a p\ue sinusoidal control motion. Consideration was given to changing 
the reduction gear ratio and thus providing more torque at low frequen- 
cies, but this was discarded because of time requirements. The gyros 
were properly iitilized, and their limitations were never exceeded; and 
the original estimates of the control deflections required to produce the 
desired responses proved accurate, making major modification of the 
original instr\xmentation unnecessary. 



Several difficvdties were encotintered during the flight test program. 
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The initial flights were conducted with the rudder locked zind with the 
ailerons driven by the oscillator. At oscillation frequencies well 
above the dutch roll the steady state condition was rapidly reached and 
satisfactory data coxild be obtained before the spiral mode was excited. 
At frequencies near the dutch roll the spiral mode was frequently ex- 
cited before the steady state condition was reached necessitating ter- 
mination of the run before any data could be obtadned. At frequencies 
below the dutch roll the initial control deflection was sufficient to ex- 
cite the spiral mode^and no data at all could be obtained. After some 
experimentation it was found that most of this difficulty could be elim- 
inated by starting the oscillations at a high freq\iency and then reducing 
the frequency to the desired value, and this procedure was adopted in 
all subsequent flight tests both for aileron and rudder oscillations. 
During the rvins involving aileron oscillations there was very little indi- 
cation to the pilots when the dutch roll frequency was being approached, 
and it was not until the data was reduced that it could be seen that the 
responses peaked at this frequency. During the runs involving rudder 
oscillations, however, the responses became quite violent at the dutch 
roll and apparently very sensitive to airspeed. During the first rudder 
run at the dutch roll frequency the airspeed control was poor, and as 
speed decreased by as little as 6 mph the responses decreased in a 
very noticeable manner. As speed was increased to the desired value 
by retrimming the elevator^ the resonant condition was again reached 
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and the responses became large. A similar variation was noted with an 
increase in airspeed over the desired value, bvit the effect was not as 
pronounced. The difference between the effects was explained by refer- 
ence to the final response cxn*ves which show a rather flat top at the 
diatch roll frequency. A decrease in velocity corresponds to an increase 
in T' and this in turn corresponds to a decrease in the dutch roll fre- 
quency. Since the actual frequency of the oscillations was maintained 
constant, the operating point at the decreased velocity was actually 
above the dutch roll. Since the experimentally determined response 
curves feill off more rapidly above the dutch roll than below it, the con- 
dition was satisfactorily explained. Because of the violence of the 
responses near the dutch roll and their sensitivity to airspeed great 
care was necessciry to keep the speed at its desired value^and rims 
where speed control was poor were discarded. 

Ref. 1 discusses some difficulty in the design of satisfactory con- 
trol locking and trim devices for the oscillation mechanism. The orig- 
inal design accomplished here utilizing the recommendations of Ref. 1 
proved highly successfxil and no diffictilty was encountered either in 
obtaining proper trim or in releasing the oscillator or the locked control 
surface when the attitude of the aircraft made this necessary. 

While it was felt in the original design that a measure of the con- 
trol field voltage to the amplidyne would serve as a satisfactory indica- 
tion of the frequency of oscillation, it was found that this assumption 



34 



was not wholly justified, A rough determination of frequency co\ild be 
made from this voltage reading, but because of residual magnetism in 
the control field of the amplidyne, the desired accviracy could not be 
attained. Precise frequency determination for data analysis was made 
from the recorded traces, but since a fairly accvirate knowledge of fre- 
quency in the cockpit during the flight testing was desirable it is recom- 
mended that in subsequent tests of this nature a tachometer generator 
be incorporated on the driving motor. 

In the analysis of the data at low frequencies it was found that, 
although the inputs were not purely sinusoidal, the steady state condition 
was rapidly reached. Harmonic analysis of both the responses and the 
forcing function was considered but was discarded as being too time con- 
suming, Since all the data had been recorded on magnetic tape, a tech- 
nique of filtering both the forcing function and the aircraft responses 
with a low frequency band-pass filter was developed. The procedure, as 
previously described, was really nothing more thAn a fast and simple 
method of performing harmonic analysis of the traces and required that 
the data be recorded on magnetic tape. The technique was a valuable 
one, and it was felt by the authors that it could be utilized . to advantage 
in all forms of steady state dynamic flight testing where telemetering 
and magnetic recording of data is performed. 

The xise of telemetering equipment in a flight test program of this 
sort requires more complicated instrimientation than the simple SFIM 
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recorder utilized in Ref. 1, in that a ground station cind its associated 
eqviipment is reqviired. The actiial instrumentation in the aircraft is 
little complicated by the substitution of the telemeter transmitter for 
the in-flight recorder. The utility derived from having magnetic tape 
recordings when the telemeter scheme is used is enormous, and it is 
felt that the advantages of this method outweigh the complications pre- 
sented by the extra instrumentation required. The filter technique of 
covirse would be impossible without tape-recorded data, but in addition 
tape records of this sort may be expanded to any desired time scale 
making phase angle measurements easier and more accurate. It is 
also possible to amplify small responses to distingviishable values, but 
care must be tadcen not to exceed the accuracy of the record. The tele- 
metering scheme vitilized here was accvirate to 1% of full scale and thus 
small responses could be expected to have a much greater percentage 
error than large ones. It was therefore desirable during this investiga- 
tion to obtain responses which were as large as possible without exceed- 
ing the range of linearity of the equations of motion. This was accom- 
plished by varying the forcing function amplitudes over the frequency 
range to give desirable responses. In general, large control deflections 
were required at high frequencies and small ones at and below the dutch 
roll frequency. 

The method used in this investigation to determine stability deriva- 
tives from frequency response was basically sound. As mentioned 
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previously, considerable difficulty was encountered in the determination 
of the cross derivatives, and it was felt that accurate determination of 
these would require much more accurate frequency response data than 
could be obtained with the simple instrumentation utilized here. The 
accurate determination of phase angles from data obtained using the 
forced sinusoidal oscillation technique is difficult indeed, eund phase angle 
errors can have a large effect in the equations for the solution of the 
stability derivatives, especially when the coefficients of these equations 
involve sines of angles near 0° or 180° or cosines of angles near 90° 
or Z70°. The technique developed here, which involved fairing a 
smooth curve through the experimental points and taking data for the 
computation of the derivatives from the smooth curve, eliminated some 
of this difficiilty, edthough it was felt that the method of averages could 
probably eliminate most of this random error even if individual data 
points were used. 

The equations derived from the rudder data were, as previously 
mentioned, ill-conditioned for solution, and this fact reqixired deternain- 
ation of some of the derivatives from the aileron data before the rudder 
equations could be used. In theory, it is possible to determine all the 
derivatives, except of course the aileron control derivatives, from 
rudder data alone. While this was not possible for the Navion aircraft 
at the condition investigated, it was felt that for other aircraft or other 
flight conditions the difficulties encovmtered here might not be present. 
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It became apparent to the aiithors that a need exists for extensive study 
of the determination lateral stability derivatives from frequency re- 
sponse data, especially as concerns ill-conditioned equations and possi- 
ble methods for their solution. While it finally proved possible to 
obtaun fairly accurate derivatives in this investigation, certain devia- 
tions from the basic method were required. The methods applied here 
were pecxiliar to this particular aircrait and flight condition and might 
very well prove unsuccessfvil for other aircraft or flight conditions. The 
best procedure to follow in the extraction of stability derivatives from 
frequency response data is, at the present time, a fxmction of the type 
and accuracy of the data, although the authors feel that with time a 
standardized procedure for eliminating difficulties could be developed. 

The use of the method of averages for eliminating random error 
instead of the method of least squares, as outlined in Ref. 4, proved 
quite satisfactory and had the added advantage of reducing the computa- 
tion time considerably. When using this method care must be taken 
when grouping the points for averaging. If automatic computing facili- 
ties are available and large amounts of data require reduction, the 
method of least squares would probably prove more satisfactory and 
less subject to error. 

It was difficvilt in this analysis to acqxiire a feel for the accuracy 
of many of the individual derivatives. As previously stated, the nature 
of the equations and the fact that small differences of large quantities 
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were involved in the solution for and S made the 

accviracy of these derivatives questionable. In the solution for , 

Ci^ , and the controlling terms involved Cy^ , amd 

CxXp respectively, and thus these control derivatives were ailso of 
doubtful accuracy. Lateral acceleration data would have eliminated 
this difficulty in the case of the side force derivatives, but was not 
obtained here. In the determination of the other derivatives, however, 
the potential degree of accuracy was not evident from the equations. 

A solution of the chairacteristic eqviation using the experimentally 
derived derivatives resulted in a period of 2. 38 seconds and a time to 
half amplitude of 1.69 seconds for the dutch roll mode. Trcinsient data 
obtained for the seime aircraft at the same flight conditions indicated a 
period of 2.40 seconds and a time to half amplitude of 1.70 seconds. 

In view of this close agreement, it was concluded that Cn^ . which 
primarily determines the period of the dutch roll mode, omd Cnj,. the 
primary damping derivative, were accurately determined. 

A compairison of the frequency responses derived from the experi- 
mentally determined derivatives and the experimental points indicated 
certain regions of discrepancy. In Fig. 19 the calculated amplitude 
response in roll to aileron oscillations was noticeably low at and above 
the dutch roll frequency, aind the calculated phase amgle was slightly 
high in the same region. An analysis of the transfer fimction contained 
in Table V indicated that am increase in the coefficient of the first power 
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term in the numerator would remedy this condition. The most obvious 
means of accomplishing this was an increase in the value of • 

An increase in this derivative would also increase the constant term in 
the nxxmerator of the yaw transfer function and thus increase the ampli- 
tude response at low frequencies. Fig. 20 shows the need for this cor- 
rection. While a detailed analysis of the effect of increasing 
was not performed, these considerations led to the conclusion that the 
value determined for was low, cuid a comparison with the results 

of Ref. 5 indicated this to be the case. 

In Fig. 22 the calcxilated amplitude response in roll to rudder 
oscillations is substantially less than the experimental response at low 
frequencies and the phase angle is consistently larger for cill frequencies. 
An analysis of the transfer function for this response indicated that an 
increase in would tend to correct both these discrepancies. In 

addition, the calctilated low frequency amplitude response and high fre- 
quency phase response for sideslip, as shown in Fig. 24, would also be 
altered favorably by a larger . A comparison with the results 

of Ref. 5 indicated that the value determined here for Cn was too 
low. 

Since the match between the calculated euid actual responses was 
generally good with the exception of the cases mentioned above, it was 
felt that the determination of those derivatives with the greatest effect 

on the frequency response, i.e., Cj • ^ , and , was 

/9 S "O ^ 
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good. While detailed curve fitting of individual responses by altering 
individual derivatives was not considered within the scope of this inves- 
tigation, such a technique for correcting derivatives might warrant 
further study. 
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CONCLUSIONS AND RECOMMENDATIONS 



The method utilized in this investigation to determine the lateral 
frequency response of an airplane proved entirely feasible and could 
xmdoubtedly be applied to many current production eur craft. The 
method, however, had the disadvantage of requiring more flight time 
than currently popular transient response methods, and thus would not 
be practical for missile testing. 

Accurate determination of frequency response previously required 
an elaborate system for the generation of an exact sinusoidal forcing 
fvmction or time constiming harmonic analysis of the data. The filtering 
technique developed here eliminated these difficulties, and could have 
wide application in this type of flight testing. 

The method utilized here for the determination of stability deriva- 
tives from frequency response data proved accurate for those derivatives 
with significant effect on the frequency response. It was, however, im- 
possible to determine Ctip . ^Ij. « f and Cl^^ with a high 

degree of accviracy. Better results could have been realized in the 
determination of Cy^ and if lateral acceleration data had been 

obtained, and it is recommended that such measurements be made in 
future tests. It is further recommended that a study be made to develop 
a systematic method for eliminating the difficulties encountered here 
with ill-conditioned equations. 
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TABI£ I 



PHYSICAL CHARACTERISTICS CF NAVION 



A. WING 



1. 


Total Wing Area, Sy, (includes flaps, ailerons 
and 19o87 Ft3* covered by the fuselage) 


l8Uo337 Ft.2 


2o 


Span bv 


33o378 Ft. 


3. 


MAC 


68o352 in. 




Angle of Incidence 
Root, ip 


♦20 




Tip, it 


-10 


5. 


Twist 

Aerodynamic 
Qe one trie 


20 31 » 

3° 


6. 


Airfoil Section 
Root 
Tip 


NACA UUl^ 
NACA 6U1CR 


7o 


Aspect Ratio, AR^ 


6.0UU 


8o 


Taper Ratio, X ^ 


0.5265 


9. 


Dihedral 


7o5® 


10. 


Root Chord 


7o2 Ft. 


Uo 


Tip Chord 


3o92 Ft. 


AIIBICN (one Aileron only) 


1200-1200 


(Beaded ekin, Frlse Ncee Balance, 
No Trailing Edge Extension) 


Wheel Throw 
Fixed bend tab 
right aileron. 


1. 


Area, Sa 


2ol61 Ft.2 


2o 


Span, b^ 


61.987 In. 


3. 


Allarcn Deflection, ^ ^ 


30° DP, 200 DN 


Uo 


Boost 


None 


5. 


Trim Tab 


Fixed bend tab 
right aileron. 



B. HCRIZONTAL TAIL 
1* Total Area 

(Includes 2.368 Ft.2 covered by fuselage) 
2 • Span f bg 
3e Airfoil Section 
U. Taper Ratio, 

5o MAC, Ct 
6. Aspect Ratio, ^ 

B-lo HORIZONTAL STABILIZER 
1. Area, Sg 

(with reference to FRP) 

B-2. E1£VATCRS 

(No Trailing Edge Extensions) 

le Total Area, Se 

2o Span, bg 

3. Deflection, ^ g 

lio Trim Tabs (32* Span, h-l/2" Chord) 

5o Root Chord 
6. Tip Chord 

C. VERTICAL TAIL 

1. Total Ax^a, Sy 

(Includes 2 .$7 7 Ft. 2 blanketed by 
fuselage and excluding 1.U83 Ft.^ 
cf dorsal fin) 

2o Airfoil Section 
Root 
Tip 



1*3.051 Ft. 2 

13.172 Ft, 

NACA 0012 
.67 

3.3li Ft. 

1*.02 

28.953 Ft.2 
- 1*0 

Smooth Skin 
Flat sided 
No trim bungee 
Balance spring 

11».098 Ft.2 

73.582 in. 

300 DP, 20° DN. 

30° UP, 300 DN. 

1.5 Ft. 

1.0 Ft. 

12.925 Ft.2 



NACA 0013.2 Mod. 
NACA 0012-61* " 



C-1. 


VERTICAL STABILIZER 






1. Area, Sv 

(Includes 0*Ui27 Ft. 2 blanketed by fuselage) 


6.873 Ft. 2 




2.C^(Wlth reference to fuselage 0 ) 


20 Nose left 


C-2. 


RUDDER 


Smooth Skin 
Rigged 3® Rt. to 
Fin 0, Fixed 
Bend*”Tab. 




AVBAf Sp 


6.052 Ft.2 




2. Rudder Deflection, ^ j. 


17® L - 23® R 




3o Rudder Pedal Throw 


5o75 in. 




Uo Trim Tab 


Fixed Bend Tab 


D. 


mOPELLER 


Hartsell 




lo Activity Factor, A.F. 


100 




2. Diameter, d 


86" 




3. Pitch G @ .75 R 


210 30 « 


Eo 


NOSE GEAR 






1. Travel 


+ 300 




2o Strut Angle (To Vertical) 


10® 


F. 


HISCELLINEOIIS 






lo Length, overall 


27o 25 Ft. 




2o Tail length, ^vT 


16.88 Ft. 


0. 


POWER 





1. Engine 



Continental E-185 Engine 
185 HP 0 2300 RPM, 29" Hg 
M.P. 0 Sea Level 



TABUS n 



EXPERIMENTAL FHEQUENCI RESPCWSE DATA 
TO SINUSOIDAL AIIERCN OSCILUTIONS 



T.A.S. 

Density Alto * 
Gross Vrto «= 



129o 6 m.p.h. 
^000 ft. 
2870 lbs. 




TABIE in 



EXPERIMENTAL FREQUENCY RESPONSE DATA 
TO SINUSOIDAL RDDIMl OSCILLATIONS 



T.A.S. • 129.6 m.poh. 

Density Alt. - 5000 ft. 
Gross Wto “ 2870 lbs. 



60 


1$ 1 




111 








1.05 


1.222 


100 o2 


7.UI 1 


176.2 




1.175 


6.1 


1o315 


1.121 


9b.5 i 


o90b 


168.li 


1.130 

‘ 1 


b.b 

1 


1.615 


1.100 


9i».6 


I0O53 ’ 


17I08 


1 1.175 ^ 


6.1 


1.85 


1.106 


79.I4 


lo2U 


■ 160.0 


1.51 ^ 


3b8o9 


1.87 


lo352 


82.5 i 


lo3bl 


I 16U.0 


1.585 


3b7.9 


2o21 


lo358 


57ol 


1.810 


1 m.o 


: 2«^35 1 


1 331.b 


2oUl 


I0I465 


37.3 


2.165 


} ,119.9 


' 2.55 

I j 


308.1 


2.U8 


I0I48 


j 5I406 I 

1 1 


1 2.13 

f 


I 132.5 


i 2.60 i 

1 1 


1 313.9 


2.58 


lo55 


1 2I06 ■ 


! 2.21 


; 106. U 


i 2 o63 ! 

< 1 


295.5 


2o58 


1.51 


1 35ob 


i 2o368 


j 12b. 0 


' 2.b6 1 


308.6 


2.69 


1.5U7 


1 -5.2 


2.218 


1 79«5 


i 2.695 j 


281.b 


2.99 


loOO 


i -2bo5 


; 1.661 


: 73.2 


' 1.99 1 

/ 


2b6.2 


3.71 


.152 


1 -51.8 


[ .822 


i 38.9 

J 

29.2 


: .980 ! 


22b. 9 


UoOO 


.29U 


! -63.6 


i .599 


.655 ' 


210.3 


i U.35 

1 


.253 


i -58.7 


! .b93 


20.0 


.552 

1 * 


217.6 


i U068 


.181 


; -66.6 


.390 


30.9 


.390 1 


226.5 


1 l*.70 

i 


! .1614 


; -66o5 


. .36b 


bl.5 


.b08 


215.7 


1 i».83 


.160 


1 

1 -68.0 

1 


.338 


20.0 


i 


- 


I 5o3l4 

i 


.lOb 


1 -77.6 


.253 


' 19o6 


.27b I 

J 


215.2 


i 5e38 

1 


* 


‘ 


; .2b35 


22.6 


.2bl 1 


215.8 


5.71 


j .083 


I -75o5 


.216 


■ 3bo5 


.2b3 ^ 

i } 


208.7 


7.21 

1 


; .0U47 


-83 o5 


; .127 


2b. 0 


^ .I5b5 i 


197.5 



TABI£ IV 



INITIAL THEOIETICAL AND EXPERIMENTALLT DERIVED 
UTHtAL STABILITT DERIVATIVES FCE THE NAVICN 



T.A.S. ■ 129o. 6 m.p.h. 

Density Alt. - $000 ft. 

Gross Wt. » 2870 lbs. 



r 

Derivative 


Initial 

Theoretical 


Results 
of Ref. 5 


Final 

EKperimsntal 

Results 




-0.592 


-0.592 


-0.591 




+0.102 


+0.131 


+0.119 




-0.0811* 


-0.0595 


-0.0513 




-0.1*60 


-0.1*30 


-0.1*31* 


""r 


-0.0937 


-0.0951 


-0.0908 




+0.085 


+0.101 


+0.0895 


S 


-0.051 


-0.0705 


-0.106 




+0.123 


+0.11*1* 


+0.206 




-0.101 


-0.137 


-0.108 




+ 0.0152 


+0.0069 


+0.0128 




+0.0138 


+0.0031 


+0,00376 

i 




-0.0625 


-0.0771* 


-0.0673 




TABLE V 



TRANSFSl FUNCTIONS FOR THE NAVION DETERMINED 
FROM EXPERIMENTAL FREQUENCY RESPONSE DATA 

T.A.S. ■ 129.6 ra.p.h. 

Density Alt. ■ 5000 ft. 

Gross Wto ■ 2870 lbs. 

Q _ -.2^^ - Z0,0^\ 5 - I 5 5* 

5 -4- <=I.OOZ 5^ + 15 . S' 4,6 5^4- 57. 4>3l 




-23.470 5^ - /7. 2-brS - /St.S'?? 

S4 + 9 . 00 ZS- 4 I3.56«^S'* 57.S-M ‘S-/.&3I 



-yp __ ,^44S^4- /I ■ ^ 3 8 S ^ 4 3 .9 86 S - ZZ.S’‘i4 

5 ( 5 ^+ ^.OCZ5^ 4 /S.Sfe? 5 ^ + 5-7. 5 -^43/) 



__ . 565 'S^ 4 - / 1 . 94 .^ 5 ^ 4 ~/.450 

5 ^ 4 9 C 023 ^ 4 + S'l.S'S'-^ S - t,Li{ 



0 

cir 



2.7935^ - 17.475" 5 - 2P.37Q 

5^ 4 9.4>oZS 3 + I3.54f5^45 7.f'f45-/-43l 



2^ — -4-376> 5 ^ - 33. 749 5*^ 4 3S-79 ^ ^ ~ 

(iv- 5 ('^T ^4 4 3 ^ 4 3 - 7 .<P(fV 5 - 






'i 








FIGURE 2 - TEST AIRCRAFT 
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FIGURE 7 - MECHANISM FOR THE SINUSOIDAL OSCILLATION OF AILERONS OR RUDDER 

WITH PROVISIONS FOR THE LOCKING OF THE UNDRIVEN CONTROL 
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Figure 1 1 - Samples of Paw Data 
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Figure 12- Samples of Filtered Data 
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CALCULATED FREQUENCY RESPONSE IN ROLL 
FROM EXPERIMENTALLY DERIVED DERIVATIVES 
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APPENDIX A 



Development of Vector Equations of Lateral Motion 
for Determination of Stability Derivatives 

The three linear equations of Lateral Motion may be written in 
the following form: 

Side force: (1) (l^ + } /^ " 



RoU: 


(2) 




c/^T) ) ■?> - 1 


/Ip pi 

" Ik 




c4^.^ 


Yaw: 


(3) 




CnpD) (p + 




C^^D) "V" + 






The coefficients of the 


se equations are defined as follows: 




C * 


- 


^yg(~zT) 




■ '^y/r 


Gi-j 




Cl* 


- 


cl(^t~)^ '3/v 




■ *=l<(r 


G-jtO 






m 












%* 


cz 






■ 






Clr* 


C3 


i zJkT ) 


r * 


“ ^n r 












V / 




r * 




‘S ( jf’r’ j 












a 


^ ( 2 .J 2 r) 











I 
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The vector technique is applied to the equations in the following 
manner. If an input to a linear system is considered to have a 

sinusoidal variation of frequency 60 , it can be considered to be a 

vector of unit magnitude lying along the real axis of the complex plane 
as shown in the sketch below. 




The first derivative of the sinusoidal input is obtained by multi- 
plying amplitude by 6J and rotating the vector through 90°. The 
steady state response of the system may be considered as a vector of 
amplitude ratio R and having a direction displaced from the real axis 
by a phase angle, . Thus, in the complex plane, the vector for 

sideslip may represented as 
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In complex notation the derivative of either the control input or the 
response for the steady state condition may be determined by multiply- 
ing by ^ oO and the integral by dividing by . 

The following are defined; 



- 


\t\ ^ 






1 # |s.n d e,^ 


A'xp ■= 






B ^ 


1 "f~ 1 

* Of 1 


A (p - 


t ar * 




B ^ = 


1 



In this notation with the input vector having a magnitude of unity the 
responses and their derivatives may be defined as follows: 



(4) 


4 - 




(5) 


D.« ^ 


- - <A> -c lo A (3 


(6) 


T(S 


= —LV "^A/3 - 


(7) 




~ A Cp T 


(8) 


D0 


~ *■' tA) B f A. 'OJ A 


(9) 


D"(P 


— -UJ ^ A0 ~ c 6L) ^ B <|) 


m 




== A-if +• 


(11) 




^ + -a.ujA''^ 



- uj A -y:- - ^ B -f- 



( 12 ) 






fl 



■* 
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I 



I 



I 
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Substitution of eqxiations (4) to (12) in equations (1), (2), and (3) and 
separation into real and imaginary parts yields the following six 



simplified equations: 

Sideforce: (13) Cy 6 Bg — Cu Bqj 



Roll: 



Yaw: 



(14) Cy'fl Aa - c* Acp = -h ou (?J>^ j 

(15) Q/ Be Cel<A<s;> ^ Ays = 

(16) A/s - uj2-A> -f 6 u3> - 

(17) -Cy,* Bs - u. A4 + ^ A ^ 

(18) -C>,^ A^ 4 ^np lO -<^nT' ^0 -4r ^ 



i 2r 



In these equations the A’s and B's will be available from the 
frequency response curves for any chosen value of to , In theory it 
is seen that equation (13) may be solved for from the response 

at a single frequency since is known. Eq\iation (14) may then 

be solved for using the computed value Cy^ . The complete 

solution of the yaw and roll equations may be accomplished if the 
responses are known at three different frequencies since the maximum 
number of iinknown coefficients in any one equation never exceeds 
three. The method of averages attempts to eliminate random error in 
single points by writing each equation for a number of frequencies and 
grouping the equations thus determined in as many groups as there are 
\inknowns. The coefficients for each group ajre summed eind the result- 
ing simultaneous eqxiations are solved for the unknowns. 

In the roll and yaw eqxiations it is found that better solutions are 
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obtained if the grouping of points is performed in such a manner that 
all the points in a given group lie in the same frequency range and 
fiirther that the three groxips chosen include data below, at, and above 
the dutch roll frequency. 

The final computational equations then become; 



Ujf\ 

Side force: (19) ^ '6^ ~ ^(A^+Au) 



(20) Cyfi ^ A/«, - Cif ■*" ^ 



uJt _ ^ '<•1 

Roll: (21a) ^ ■+ ^ u; 









IT V At 

^ Txi .I'K, -Cj^ <toA> 

(21b) 4- Qp £tuA(p+ — 



tVjL 






(21c) 






^ (U.V , -p 1 

(22) fAA - ^ ~~ f u.M-4 + 






Yaw: (23a) - C./ £ Bfl " 1 



tJU (L 



cot 






U.^-vrw 



a.^.v*v -r Ai — 'Cuj^R^/ 

(23b) - r,;: f S6 - £- . c.., - 2^ H 

iug, 



X>L? <io^b<() + <^>vr fi/ 



(23c) - ^n,s ^'nf, £wA4)+ 

^ . U.).^ — ;,. ... ^ / 9 

- 2A« + Cp *■ fl f,w‘A<{. -Cnr 2:“'AV--'!‘,.^= ^‘■'A;: 






U>^ 









( 24 ) 
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The development of the computational equations was performed 
for the rudder case with the coefficients Ci <■ and C_ , equal to 

“ a<rv- 

zero. It is obvious that for the aileron case the equations are identical 

to those developed above, with the substitution of and ^ 

for Cl f and Cn r and C„ ; 0 . 

Ar y 

If difficulty in the solution of the equations ajrises because of ill- 
conditioning, it may be possible to resolve in some direction other than 
that of the control deflection vector. It may also be possible to elimin- 
ate some troublesome coefficient between equations and thus obtain 
solutions. No attempt was made here to consider all the possibilities 
of modification of the equations. 

If the derivatives are known, the transfer function coefficients 
may be easily developed as follows: 
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- C* 5 






*■ 


sVG^S 

lx 




* 

j 

i 

J. _ J 


1 (2 * 
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* A /• C ^ 

- S ^ ^ r ^ 




Sr 


S T Cy C 


c 5 

- u 












1 


\ ~ 


JX 2 .. 2 . p X X n *'<. 

— - J "" ^/v\p s + ^ r ^ 


6 _ 


-t- 


C^S'^'fC^S + Cg 





+ Ci. 5 -e C. 
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5 4- Cn/ a 


c ^ 


5 




* 

C Xfi 




II* 
Ix '' 




~ Ctx/6 


C.-J* 


S*4 Cr>r 5 






- cv* 


S 






5 02^ ^ 


■^ 1 5 ^-f- Cil^ 5 

Xx 




- C«,3* 


— S 


5^4 Cv.7 5 


^ zi 


C 9 5-2 +. 


C 1 0 ^ ^ \ \ 




<^r 


5"^ 4- c, 53 


+ Ci5"-f C 3 S 


■+• c^ 




4^• 

5 4- Cy ft 


Cl* 


cy^:* 






5^4 C 2 * 5 






-Cnfi* 


-Xx?,2. f- * -r, 
. j - J 

4? 


Tt 

C^Jr 



C+ Cv/^ 


Cl.* 


5 


Cie* 


5^4 CH^S 


^^»<=>- + C£p5 
Tx 


~ C 


~ => ^|S“ 

X? 


5^4 C^* S 



' 4 ' 
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C/2. 5 C/ 3 5 0-4 






5 (s^ C, 5 ^ -h Cl 5 ■* 
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The coefficients of the transfer functions are determined as 
follows, where Cy^ ~ • 



C. =: 



ix It 



c,' = Cy* + c4 ^ C 






* hi. - CJ* 



Ye 



* r. ^ 



c, = Cn,: <• c/^ a; + cyj c„* + q; - o; - c.; Cn; 



T*- „ y(- 



* r. 



I'x r.v: 

^ 4" 






ix a 



~v>_p 



^ •4 

- c -''' 



Ix B 
Yft '--^r x7 



C 3 - Cl 4- Cxp C-^ - 0_4> Cvfp + Cy^ Ci.p Cn,- - Cy^ C 






Ikb p ’*’ 

Cx| = Ci^ Q/3 C.Kr "" Cl Cp.^ Cr.fi 



Cs 



/ 

^'/ir { 






Cvi*Cl^ + Cy;,C*y - 1|-| Oj" - C„'^ - 



- Cy/, CV Iff 



c,‘ = Cy/, c4c,%. - c«*, c4 + Cv>^,CiT 

+ “ ^Y^r 

cJ = + a/, c„;c; 



C<5 = + Cp y 






i IT T 



. X. 



c 



(O 



Cy^ + Cy^’^^ Ll? - Cy/^C^’^ -h 

+ Cn^t- Cy^ 



•3f 

' r 
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C, 



3 Cv\^ + Cv(^ 4 ^^'nf C -^Q + ’^'^if^y\.$ 



4 If T 

~ ^ f. a C-ri f. + C'/ /S 



C.,. 

Ci5 



C 



^ >f 

V' 



•Tt 

■^Sr 



Cy,\ + £^| CJlT^ 



- Cn^^C^(3 +- Cv6 - 04 



4 Cy. C 



£ 



X* 






•f Cnp 4 *^Yfi 



I/? 

T£ 



C,4 = 



- Tt 



15 



‘'V^ C^.p - ^'{^f04is Chp + Cy^^Cip Qk/S 
-f Cya C-»xp 



= cmIcIc;!:^ 4- Cntc*c4 



Once the transfer fimction coefficients are known the frequency- 
response maybe determined by the substitution of Llo for s, giving 
the following: 









e4CfcLo^)^-h (c^uj - 

- C C 4 ) 4- (CjU^ 






t,,-' 1 _ /C,w - C.U.3 
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CO 
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1 iv- 





(C(5 (^(4^ -C,2,Uj~) 

(C j to ^ - C 5 to^ ) to ST _ -*- a?) 
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ilClirt 



Ci^to— Ctitu^'i 

C 15 — C Lw 2, ] 
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